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0\ '. Abstract 



^O I We study numerically the inhomogeneous pre-big-bang inflation in a spher- 

ically symmetric space-time. We find that a large initial inhomogene- 
ity suppresses the onset of the pre-big-bang inflation. We also find that 
even if the pre-big-bang infiationary stage is realized, the initial inhomo- 
^L I geneities are not homogenized. Namely, during the pre-big-bang infiation 

(^ I "hairs" (irregularities) do not fall, in sharp contrast to the usual (potential 

^^ ' energy dominated) infiation where initial inhomogeneity and anisotropy are 

f"^ ■ damped and thus the resulting universe is less sensitive to initial conditions. 

§^ ; PACS numbers: 98.80.Hw, 04.25.Dm, 04.50.+h, 98.80. Cq 
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I. INTRODUCTION 

The pre-big-bang(PBB) cosmological scenario inspired by string theory has been sug- 
gested as a possible working mechanism of inflation within the framework of the low energy 
effective action of string theory [|l|. In this scenario, the inflation occurs by the kinetic energy 
of a massless dilaton, and the curvature and the string coupling grow until the effective action 
breaks down. After that epoch stringy nonperturbative effects become important and hope- 
fully enable the universe to make a smooth transition to a standard Friedmann-Robertson- 
Walker(FRW) epoch of decelerating expansion (for recent attempts toward graceful exit, see 

One of the aims for the introduction of the inflationary universe scenario ^ was to solve 
the homogeneity problem. However, most of works on the PBB inflation is done within 
the context of a homogeneous space-time (however, see 0]). So, questions arise as to a 
homogenization process due to the PBB inflation: Does the PBB inflation really homogenize 
initial inhomogeneities? Can initial inhomogeneities prevent the onset of the PBB inflation? 
The purpose of this paper is to answer these questions in a spherically symmetric space-time. 

The peculiarity of the PBB inflation is that even classically the weak energy condition 
1^ is violated during that epoch when the equation are written in the form of the Einstein 
equations ||^, while the strong energy condition is violated but the weak energy condition 
is respected classically in the case of the usual (potential energy dominated) inflation. This 
suggests that the behaviour of inhomogeneities in the PBB inflation is quite different from 
that in the usual inflation. 

According to the linear perturbation analysis, the density perturbation 6 oi p = wp 

2(l+3u)) 

matter, with w being a constant, behaves as (in the long wavelength limit) 6 oc 1 3(1+™) . 
Similar behaviour is found by the approximation method to describe the super-horizon 
scale inhomogeneity (the gradient expansion method) p|. Therefore, for the strong energy 
violating (but the weak energy condition respecting) matter (—1 < w < —1/3), the density 
perturbation decays with time. An inflaton generically satisfies this property. This behavior 
is related to the so-called cosmic no hair conjecture, which states that all expanding universe 
models with positive cosmological constant asymptotically approach the de Sitter space-time 
0. On the other hand, however, for the weak energy condition violating matter {w < —1), 
the density perturbation grows with time. 

Of course, the dilaton does not have a constant w and so the above consideration is 
not directly related to the PBB inflationary cosmology, but it strongly suggests that the 
behavior of the density perturbation in the PBB inflation is quite different from that in 
the usual inflation. In order to study how initial inhomogeneity affects the onset and the 
duration of the PBB inflation, we need to consider large inhomogeneities, and hence we 
have to resort to fully numerical treatment. In this paper we investigate the inhomogeneous 
pre-big-bang inflation in a spherically symmetric space-time by solving the field equations 



fully numerically. A related numerical work based on a spectral method was appeared [[10 



however, the analysis is limited to weakly nonlinear perturbations. There also appeared 



a related analytical work ||rT| where a criterion for PBB inflation is proposed, which is 
consistent with our numerical results. 

This paper is organized as follows. In Sec. 2, basic equations based on the ADM formalism 
are given. In Sec. 3, after the details of the initial condition and the numerical method is 



described, the numerical results are given. Sec. 4 is devoted to summary. 

II. FIELD EQUATIONS 

The low-energy effective action is given by 

S = l d^xy^^ ( (^)i? + ^^'^ V^0V,0) , (2.f ) 

where Is is the string scale which is of order the Planck scale, is the dilaton, and ^^^i? 
is the four dimensional Ricci scalar. We have omitted other matter degrees of freedom for 
simplicity. Varying the action by g^'^ and (p yield 

G^u = -V^V,0 + ^(7^.(V0)2 = T^„ (2.2) 

□0 = (V0)^ (2.3) 

where V^ and □ are a covariant derivative and the d'Alembertian of Qf^^, respectively. 

Although after the conformal transformation such that QEfiu = g'^O/iu the action is 
reduced to that of a massless scalar field coupled minimally to the Einstein gravity, we 
shall work consistently in the string frame. This is because the gauge condition in the 
string frame metric has a direct geometrical meaning and because the interpretation of the 
numerical results is direct and clear as suggested in Wl 



A. Basic Equations in (3+1) Form 

In the ADM formalism, the line element generally takes of the form 

ds'^ = -a'^dt^ + ^i0dt + dx'){(3^dt + dx^), (2.4) 

where a, (3^ is the lapse function, the shift vector, respectively, jij is the three-metric of the 
spacelike hypersurface. A timelike vector tt,^ normal to the hypersurface is given as 

n^ = (-a,0,0,0). (2.5) 

The extrinsic curvature Kij is then defined by 

^ij = -^£nlij, (2.6) 

where £„ denotes the Lie derivative along the integral curve of n^. The field equations 
Eq.( p.2|) then break up into evolution equations 

7ii,i - £f3lij = -2aKij, (2.7) 



Kij^t - £i3Kij = -DiDja + a 
and constraint equations 



Ri, + KK,, - 2KuK', - Uj + ^%jipH - 5/) 



(2.8) 



R + K^- Ki, K'^ = 2pH, (2.9) 

D,K,^ - D,K = J,. (2.10) 

Here Di is the covariant derivative operator with respect to jij, Rij is the three-dimensional 
Ricci tensor, K is the trace of the extrinsic curvature, and 'jij^t denotes the partial derivative 
oi'jij with respect to t. pn, Ji and Sij are defined by the right-hand-side of Eg. ( pl2D (denoted 
as 1 np ) as 

PH = T^un^7f, (2.11) 

Ji = -T^,n''hi\ (2.12) 

Dij = Ifj^iytii hj , \Z.l6) 

where 

hf,u= 9tMu + nf,n^ (2.14) 

is the projection operator on these hypersurfaces. We also need to write the equation of 
motion of dilaton Eq. (|2.3| ) in the first-order form. To do so, we introduce the following 
variable likewise Kij |T^ 

U = -£r,<P = --Ut-/3'<P,i). (2.15) 

Then the dilaton equation of motion can be written in the first-order form as 

0,t - /?>,i = -«n, (2.16) 

n,t - /?^n,i = a (kVl - n^ + (D^)^ - D'^(^ - 7^^a,i0j. (2.17) 



B. Basic Equations in Spherical Symmetry 

Now we specialize the above general considerations to a spherically symmetric space- 
time. The most general line element of a spherically symmetric space-time is written in the 
form 

ds" = -{a^ - A^(3^)dt^ + 2A^(3dtdr + A^dr^ + B^^dVt^ (2.18) 

where (3 is the radial component of the shift vector. Because of the spherical symmetry, only 
two components of Kij are regarded as independent variables. We introduce Ki = K'^r- 
Because there does not appear "bare" in Eqs.( |2.2|) and (|2.3| ), it is convenient to introduce 
the radial derivative of it as another independent variable: 



$ = -^. (2.19) 



r 



Then the evolution equations in (3+1) form are written as 



A^t - pA^r = -aAKi + (3^rA, 

B^t - PB^r = -laB{K - K,) + ^B, 
2 r 



K,t - I3K, 



+ a 



1 / { B,j, A^ 
—77 Axa XX + Qci x + 4:xa x 2—^; '— 



-K^ -KKi + -K^,+KIi-W + ,, 

2 ^ 2 ^ A2 



A2 ^V 5 



2^ + 3 ) $ + 2x^,x 



1 / A 



+ a 



4 / S, 






+ KKi - K,U + — (2x^,x + $ - 2a; -^$ 



n,j - /3n,^ = aKU - a n^ 



x$^ 



A2 



a / xa . 
— 2 — '- 
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) 1: 



xB. 



xa 



$,t - /3$,, = -2a xU - 2aUx + Ar + - * 



+ 4^ + 3j$-2^$,, 

/3^ 



(2.20) 
(2.21) 



(2.22) 



(2.23) 

(2.24) 
(2.25) 



where we have set x = r"^ to include the regularity condition at the origin. There are also 
constraint equations 

AB ^ 1^ ~ x_ 



2 

A^ 
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+ 3 $ + 2X$ ; 



Ki^x - K,x + (^ + ^) (3^1 - ^) = -n,. + ]^K^^. 



(2.26) 

(2.27) 



At this stage, it may be instructive to recover the homogeneous and isotropic solutions. 
This can be done by taking a = l,/3 = 0,A = _B and neglecting the spatial dependence in 
the variables. Then we arrive at the following equations: 



A 3 



s* 



n, 



KU-U^ 



'-K^ = 2KU - W 



(2.28) 

(2.29) 
(2.30) 
(2.31) 



with n = —(p^f By rewriting as if = —K/3, Eq.( ^.28| ) is nothing but the definition of the 
Hubble parameter. Eq.( p.3lD is the Friedmann equation. Eq.( |2.29| ) and Eq.( p.30D are the 
evolution equations. Either of them is redundant. Eq. (|2.31| ) admits the following solution: 



K = -3H = ^^^n. (2.32) 

It is easily found that minus sign corresponds to the PBB branch solution, while positive 
sign corresponds to the post-big-bang branch solution. 



III. NUMERICAL RESULTS 

A. Initial Condition 

Initial data must satisfy the constraint equations Eq.( p.26| ) and Eq. (p.27|) . Because of 
spherical symmetry, we can consider, without loss of generality, the conformally flat three- 
metric: 

df = 'yijdx'dx^ = ■^{r)\dr^ + r'^dQ^). (3.1) 

As a choice of the initial time slice(or choice of the extrinsic curvature), we adopt the 
following extrinsic curvature considering the similarity with the homogeneous universe 

K^, = lK^^J + CT,, = ^(1 - ^)n7^„ (3.2) 

where aij is the traceless part of the extrinsic curvature. That is, we take aij = and 
K = {3 — v^)n/2. The former condition is for simplicity while the latter condition is chosen 
so that it corresponds to the Friedmann equation in the PBB phase in the homogeneous 
universe (see Eq.( p.32| )). 

Then the constraint equations become 

Ax^p,^^ + Qip,^ = xiIj,^^ + -^<l> + 1^$,^. - |^<l>^ (3.3) 

n,. = ^^^$n. (3.4) 



Note that because of the choice Eq.( |3.2| ) the Hamiltonian constraint Eq.( |3.3D is now in- 
dependent of Kij and 11. Given $, we solve Eq.( p.3|) iteratively to get ip by the cyclic 
reduction method for a tridiagonal matrix. However, the right-hand-side source term does 
not have an appropriate ^/'-dependence for numerical treatment: as ip increases, the source 
term also increases, and therefore the equation is numerically unstable. So, we "redefine" $ 
by $/ = ip^, and rewrite the right-hand-side of Eq.( p.3| ) in terms of $/. That is, 

4x7^,.. + Q^P,. = ^$7 + l^i + 1$/,. - ^$?. (3.5) 

Physically, fixing $/ corresponds (roughly) to fixing the proper kinetic energy density of the 
dilaton field. The boundary conditions are 

tr = (3.6) 

at the origin for regularity and 



{r{^ - 1))^, = (3.7) 

at the numerical outer boundary to guarantee the asymptotically flat Friedmann property. 
Eq.( |3.7| ) corresponds to setting the asymptotic scale factor to be unity. We then solve 
Eq.(|^ to get n. 

We consider the following Gaussian form for $/: 

We examine whether changing the scale of inhomogeneity (i.e. A) can affect the onset or the 
duration of the PBB inflation. We use the unit where Hi{r -^ oo) = —Ki{r -^ oo)/3 = 1 
initially, which determines the time scale to to the big bang singularity (in FRW universe) 
such that to = H~^ /\/?). The normalization factor 0o is arbitrary and we set it to one. 



B. Numerical Details 

We describe the numerical details for solving the field equations by the finite difference 
method. Time steps are labeled by the index n and spatial grid points are labeled by i. We 
define the time derivative operator 

/ V \n — ^tn-1 ^i — ^i ^t„ Zj^ — Zj^^ 

^ '''' ~ Atn + Atn-1 Atn At,, + At„_i At„_i ' ^ ' ' 

where At.„ = t„+i — t„ and a operator for a derivative with respect to x = r^ 

iZ,.)7 ^ ^^^^(^r+i - Zr) + ^^±LZll^zr - Zl^. (3.10) 

Then for a grid uniform in r, the evolution equations are second-order accurate in both space 
and time. We take a nonuniform grid because we need to perform numerical calculations 
with various initial conditions and therefore it is desirable to use small grid points. Typically 
we use 512 spatial grid points in this work. The time step is constrained by the Courant 
condition. For a nonuniform grid, we use 

At = eimin ( /^+^ ~ ''' ] , (3.11) 

and we take ei = 0.2 for accuracy. The Courant condition is sufficient for stability, however, 
not for accuracy. Therefore we impose the additional requirement that A, B, $ do not change 
too much (near origin) in one time step: 

f A B $ \ 
cit = .,min(^-,-,-j, (3.12) 

and choice of the value of £2 depends on the initial condition. For example, we choose 
£2 = 10-3 for A/Hr^ = 1.0. 

Regarding the boundary conditions, Eq. ( |2.27| ) requires a single boundary condition. We 



set Ki = K/3 at the origin for regularity. Evolution equations Eq. (|2.22|) -Eq. (|2.25|) require 



an outer boundary condition as well as the regularity condition at the origin of kind Eq.( |3.6| ). 
To determine the value of a variable at imax, we extrapolate linearly in x using the values of 
the variable at imax — 1 and i^ax — 2. This is because it is not obvious that the conventional 
out going- wave boundary condition should be imposed in the present context. Since we locate 
the numerical outer boundary at the point much further than the scale of inhomogeneity 
(> 10 X A) where the universe is approximately FRW, such a choice is sufficient in practice. 



C. Numerical Test 

As a test of the numerical code, we have used homogeneous initial conditions and com- 
pared the numerical solutions with the exact homogeneous solutions. The results are shown 
in Fig.l. There the exact solutions are shown as solid curves, while the numerical solutions 
at the center (r = 0) are plotted as crossed points. We find excellent agreement. In fact, 
the maximum relative error between the numerical solution and the exact solution was less 
than 10"'^ and the maximum relative error in the Hamiltonian constraint was less than 10~^ 
because the variables appear only quadratically in the constraint. 

Barrow and Kunze derived an exact solution of the collapse of a homogeneous spherical 
region of the stiff (p = p) fluid |T^ . The comparison with their solution is difficult in the 



present work because the spatial distribution of the dilaton is step-function-like and hence 
$ becomes singular near the surface (delta- function- like). Furthermore, Eq.( p.24| ) contains 
the spatial derivative of $, which develops a further ill-behaved singularity. 

D. Results of Numerical Calculations 



We choose the geodesic slicing condition, a = 1 with zero shift, (3 = 0. We use Eqs.( |2.2(] - 



2.22| ) and (|2.24| - |2.25D to solve A, B, K, 11, and $ and use the momentum constraint Eq.( p.27|) 



instead of Eq.( ^.23 ) to solve Ki. The Hamiltonian constraint equation Eq.( ^.26 ) is used as 



a check of the numerical accuracy. As a geometrically invariant diagnostics for the inhomo- 
geneity, we calculate the Bel- Robinson tensor [|15|. The Bel- Robinson tensor is defined in 
terms of the Weyl tensor as 

where C^upa is the Weyl tensor and *C^ypcr = \^p.uapC^^ p^ is its dual. The "superenergy" 
density on the three-dimensional spacelike hypersurface with n'^ being the unit normal vector 
is defined by 



W = Tp^p^nVn'n'^ = E^^E^^ + B^,B^\ (3.14) 



where 



Ep^^Cp^^pn'^n^, Bp,= *Cp^,pn''n^. (3.15) 

They can be written in the present case as 



Eij = Rij + KKij - 


- K/Kij 


2 {^'^ ~ 


~ll.A)- 


2 

- ^lijPH, 


(3.16) 


Bij = eim(iD K"^j). 










(3.17) 



We have numerically solved the field equations with various choice of A. Numerical results 
are shown in Figs.2,3,4 for A/H'^ = 1, 5, 0.2. 

For A/ H~^ = l(Fig.2), which can be regarded as a fiducial case, the evolution was ter- 
minated at t ~ 0.49 beyond which numerical calculations became inaccurate near the origin. 
The inhomogeneities in $ and 11 always grows during the evolution. The superenergy of the 
Bel- Robinson tensor shows that after the decay of the initial irregularities, new irregularities 
emerge and grow. This phenomena can be understood as follows. Initially inhomogeneous 
region expands being dragged by the cosmic expansion. However, its self gravitational en- 
ergy dominates over the background energy density soon, and then that lump decouples 
from the cosmic expansion and begin to collapse. We note that although irregularities are 
increasing, the universe nonetheless expands rapidly since the expansion rate is increasing 
(See Fig. 2(a) (b)) The situation is to be contrasted with the usual potential energy dominated 



infiation (for example, see Fig. 14 in |jT6[ and Fig.3 in [T] 

For A/H~ = 5(Fig.3), the evolution was terminated at t ^ 0.46. Further evolution made 
the numerical calculations inaccurate near the origin. We find that the behavior similar to 
A/H~^ = 1 case appears again. For A/H~^ = 0.2(Fig.4), the evolution was terminated at 
t ~ 0.12. Even if we perform numerical evolution further, the evolution was forced to stop 
soon. Presumably the singularity is reached at the origin first. Note that K and Ki are 
growing positively near origin. That means the expansion is delayed and eventually turns 
into the contraction to the singularity. Consequently, the metric (or scale factor) does not 
grow much near the origin. We do not see the growth of inhomogeneity in the Bel-Robinson 
tensor probably because there is no time for them to overcome the cosmic expansion. If we 
could perform the numerical evolution any further more accurate, we could see the same 
structure as in Fig.2 and Fig.3. 

We also calculate the e-folding number A^ of the expansion at the origin 

A^ = In I (A(r = 0)B{r = Off^^ li){r = Of\ , (3.18) 

and the results are shown in Fig. 5. For comparison we also plotted the e-folding number for 
the FRW universe. We find that for A ^ 0.3 x H^^^ the universe does not infiate sufficiently 
and/or recoUapses. Note that the e-folding number for A ^ 0.4/7^"^, where the onset of 
PBB- infiation is not prevented (in agreement with the rough estimate in |11[), is larger 
than that of the homogeneous solution. This is because initially the inhomogeneous dilaton 
is superposed on the background homogeneous solution and that inhomogeneous energy 



contributes in turn to the superinfiationary expansion. In ||TT[, a criterion for PBB infiation 
is proposed that is translated roughly as, A ^ Hj^^, which is consistent with our numerical 
results. 

Our results may be understood intuitively in the Einstein frame. In the Einstein frame, 
the PBB phase corresponds to a recollapsing universe with a massless scalar field rather 
than an expanding universe. So, inhomogeneities should grow. 

To conclude, we find that a large inhomogeneity, A < 0.3H~^, reduces the amount of the 
PBB infiation and sometimes even does not allow the development of the PBB infiation as 



in the case of the usual inflation |T^. However, even if the universe enters the PBB inflation 



stage, the initial inhomogeneity grows and is not smoothed out globally at all. 

E. Fine- Tuning Problem 

We consider the meaning of the condition for the onset of the PBB inflation(A ^ 0.4iJj~^) 
from the viewpoint of the resolution of the horizon problem |]TU . 



Let the time of the beginning of the PBB inflation be ti and the final time be tf and 
the time of the final singularity be ts- tf is supposed to be the time when the stringy 
nonperturbative effects become dominant and the low energy effective action Eq.( p.l| ) is no 
longer valid. The decrease of the comoving Hubble length, {Ha)~^, measures the amount of 
inflation. Using the PBB solution in FRW background |1|,0, we have 

Z*^/ _ ( ^s ~ w \ ^ fis — ti\^ , , 



Hitti \ts — tf J \ I 



where a is the scale factor and H^ ^ ~ t^ — tj ~ /^ is assumed. Z > e®° is required to 
solve the horizon problem |^ . On the other hand, we find that the lump of the initial size 



A > OAH^ ^ inflates. Using Eq.( |3.19| ) this inequahty can be rewritten as 



A > OAH-^ ~ 0.4(t, - ti) > e^Hs, (3.20) 

which shows that initially extremely large inhomogeneous region (in the unit of the string 
scale) is required in order to get the PBB inflation and to solve the horizon problem. This 
large number may be regarded as fine-tuning of the initial condition because Is is a natural 
length scale in string theory. 



IV. SUMMARY 

We have studied numerically the effect of the initial inhomogeneities on the onset and 
duration of the PBB inflation. We found that a large initial inhomogeneity does suppress 
the onset of the PBB inflation as in the case of the usual inflation [rH|. Further, even if the 



PBB inflation is realized, the initial inhomogeneity grows contrary to the usual inflation. 
The initial scale of inhomogeneity is required to be extremely large (in the unit of the string 
scale) to realize the PBB inflation and to solve the horizon problem. 

The fine-tuning problem of the PBB inflation was recently addressed in ||19| , ^ for ho- 
mogeneous cosmologies in the context of the horizon problem and the flatness problem: 
the pre-big-bang scenario is quite sensitive to spatial curvature and anisotropy. Combining 
these results with our results, it may be said that if the present universe would be resulted 
from the PBB inflation, the universe would have to be initially extremely homogeneous and 
isotropic and flat. 

Of course, our study does not exhaust all the possibihty; (i) we did not include other 
massless degrees of freedom (for example, the axion), (ii) the asymptotic condition was 
limited to the flat FRW universe, (iii) the space-time dimensionality was fixed to four. It 
would be interesting to relax the above situations. 

10 



Hawking and Moss once said in P that the asymptotic approach to the de Sitter state 
(in the usual inflation) is "very similar to the way that a gravitational collapse rapidly 
approaches a stationary black-hole state (outside of the black-hole horizon) which depends 
only on the mass and angular momentum but which is otherwise independent of the nature 
of the collapsing body" . By the same black-hole analogy, the PBB inflation may be similar 
to the "interior" of a black-hole which is dependent on the nature of the collapsing body 
(initial conditions). 
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FIGURE CAPTION 

Fig.l. Solutions with homogeneous initial conditions. Solid curves correspond to the exact 
homogeneous solutions, while crossed points correspond to the numerical solutions. 
We find excellent agreement. 

Fig. 2. Evolutions from inhomogeneous initial conditions up to t ~ 0.49 x H^^ for A = 
1.0 X Hj^^. All variables are normalized by Hi. Thick curves are initial values. Time 
evolution is in the direction from top to bottom for K,Ki,^,Il and the relative error, 
while bottom to top for A and B. For the superenergy of the Bel-Robinson tensor, 
first the initial lump decays, and then new lump appears near the origin. 

Fig. 3. Evolutions from inhomogeneous initial conditions up to t ~ 0.46 x H~^ for A = 
5.0 X H~^. Thick curves are initial values. 



'I 



Fig. 4. Evolutions from inhomogeneous initial conditions up to t ~ 0.12 x iJj ^ for A = 
0.20 X H[^^. Thick curves are initial values. 

Fig. 5. E-folding number at the center as a function of time. A/H~^ = 

5.0,1.0,0.5,0.4,0.3,0.2,0.1 from top to bottom. Dashed line corresponds to the ho- 
mogeneous solution. We find that for A ^ 0.3 x H~^ the universe does not infiate 
sufficiently. 
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